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Abstract 



We compute cross sections for the Drell-Yan process in nuclear colli- 
sions at next-to-leading order (NLO) in as- The effects of shadowing 
on the normalization and on the mass and rapidity dependence of these 
cross sections are presented. An estimate of higher order corrections is 
obtained from next-to-next-to-leading order (NNLO) calculation of the 
rapidity-integrated mass distribution. Variations in these predictions 
resulting from choices of parton distribution sets are discussed. Nu- 
merical results for mass distributions at NLO are presented for RHIC 
and LHC energies, using appropriate rapidity intervals. The shadowing 
factors in the dilepton mass range 2 < M < 10 GeV are predicted to 
be substantial , typically 0.5 - 0.7 at LHC, 0.7 - 0.9 at RHIC, and ap- 
proximately independent of the choice of parton distribution sets and 
the order of calculation. 



Introduction 



In a previous study flj we provided a systematic survey of theoretical predictions 
for the Drell-Yan process P, ^ in nucleon-nucleon colhsions at energies relevant to 
ion-ion experiments at RHIC and LHC In this study we extend our work to nuclear 
collisions at the same energies. A short theoretical discussion can be found in and 
a more complete review of the perturbative QCD calculations in van Neerven's article 
p. As before, we explore the uncertainties in the dependence of the production rate 
on the dilepton's mass M and rapidity y due to the choice of parton distribution 
set. In addition, we investigate thoroughly the effects of nuclear shadowing on the 
normalization and on the mass and rapidity dependence of these cross sections. 

Our predictions for da^"^ /dMdy are based on a perturbative analysis of the un- 
derlying partonic processes to order H, i, ^ H- The parton cross sections are 
folded with parton distributions which are extracted from data on the basis of next- 
to-leading order analysis. Unfortunately, only leading order determination of nuclear 
shadowing effects is available presently, partly due to the limited amount of data 
available for such a study. Even though it deserves a careful study, we feel that the 
main features of shadowing will not change qualitatively when going from leading 
order to next-to-leading order analysis. 

As in earlier studies, we find that the perturbative corrections and the uncertain- 
ties due to the choice of parton set grow as M decreases. From the point of view of 
the heavy ion physics, the mass region from 2 to 10 GeV is of most interest. In order 
to estimate the applicability of perturbative calculations at these relatively low mass 
values we recall results for the rapidity-integrated dcr/dM, for which the complete 
NNLO calculations (order a1) are available. The relative magnitude of the NNLO 
correction sets one limit on our confidence in the applicability of perturbation theory. 
We show in Figure ^ the mass dependence of the i^-factors from a calculation with 
NLO parton distributions and parton level cross sections up to NNLO. One sees that 
for M > 2 GeV at RHIC energy and for M > 3.5 GeV at LHC energy the NNLO 
contribution is at most about a 5 % correction to the NLO result. For lower masses at 
LHC energy the NNLO contribution grows sharply. At M = 2 GeV, the NNLO cor- 
rection alone is about 25 %, which is larger than the NLO correction itself of about 20 
%. Hence in this mass and energy range we must consider the calculations as merely 
a general indication of magnitudes, since convergence of the QCD perturbative series 
is not yet evident. Keeping this caveat in mind, we restrict ourselves in the following 
calculations to the NLO results. 

We will not repeat here the discussion in on the dependence of the cross sections 
on the renormalization scale /i^ and the factorization scale fip, or choice of MS or DIS 
regularization scheme. Since the determination of parton distributions has improved 
considerably during last few years, we calculate the NLO cross sections utilizing the 



recently-extracted CTEQ 0, GRV [jlO|, and MRS sets. In this connection we 
discuss also the dependence of the expected shadowing effects on the choice of parton 
distribution set. 

In the next section we will define the average nucleon-nucleon cross section in 
a nucleus-nucleus collision and, using the nuclear transverse overlap function, the 
number distributions normalized to one collision. We will also briefly introduce the 
nuclear shadowing modifications to the parton distributions. In the following section 
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Figure 1: K-factors at RHIC energy (200 GeV) and LHC energy (5.5 TeV) for NLO 
(dashed lines) and NNLO (solid lines) Drell-Yan calculations of da/dM using MRSG 
structure functions in proton-proton collisions. 

we present the NLO results at RHIC and LHC energies. The effects of shadowing 
are then implemented, followed by numerical results on the rapidity-integrated cross 
sections. A discussion of general properties of the shadowing correction completes 
this work. 

Cross sections and shadowing modifications 

Let us consider a collision of two nuclei with mass numbers Ai and A2. Since the 
Drell-Yan process is sensitive to the isospin structure of the nuclei, we must specify 
also the proton numbers Zi and Z2, and neutron numbers A^i and N2. We will assume 
that the shapes of the proton and neutron density distributions within a nucleus are 
the same: 



where p^(r) is the total nucleon number density. In the numerical calculations we 
use the Wood-Saxon parametrization 



PpW = Pp(^>s) = -p^(r). 





(2) 



where no = 0.17 fm~^ d = 0.54 fm and = l.UA^^^ - 0.86^-^/3 ^^^^^ 

the nuclear thickness function T^ {s) and the overlap function T^^^_^(b) at impact 
parameter b as 

^A,A.(b) = /dsT^^(s)T^^(b-s). (3) 

In a p — A collision we expect the total cross section of a hard process to be the sum 
of individual nucleon-nucleon cross sections, 

dapA = Zdapp + Ndap^. (4) 

The corresponding cross section daAiA2 for + ^2 collision can be written as 

do-^iAa = Z1Z2 do-pp + N1Z2 dcTnp + Z1N2 do-pn + A^iA^2 dcTnn • (5) 

The average number distribution in one collision at impact parameter b is obtained by 
multiplying the average nucleon-nucleon cross section in Ai + A2 collision, dcr^^^^^ = 
d(TAiA2Mi^2, with the overlap function TAiA2(b). In particular, for the mass distri- 
bution per unit rapidity in central, zero impact parameter collision we have 

dydM A1A2 dydM ^ ' 

In the QCD-based parton model the Drell-Yan cross section for nucleon(rii)- 
nucleon(n2) collision (n^ = p, n) can be expressed as 

da-.,., =^/r(^i,/i')®/r(^2,/i')®da^^. (7) 

in terms of partonic cross sections (or coefficient functions) da^j for partons i and j 
and parton distribution functions fj^^ {x , fi"^) of parton i in nucleon rii. In numerical 
calculations it is practical to combine eqs. (||) and (0) to join the parton distributions 
in proton and neutron using appropriate weights to a parton distribution in an average 
nucleon. 

The distribution of a parton flavour z in a bound proton of a nucleus A can be writ- 
ten as fi/^i^y Q^) = Q^)fi{^: Q"^)) where /f (x, Q"^) is the corresponding parton 
distribution of the free proton, and Rf{x, Q"^) defines the nuclear effects (shadowing) 
at each x and Q^. The parton distributions of bound neutrons can then be obtained 
by approximating /"(^)/^ ^ fd(d)/A fd{d)/A ~ fu{u)/A (exact for isoscalar nuclei) 
and /jy^ = for the isospin symmetric distributions. 

In practice, we include the nuclear effects in the computation of the Drell-Yan 
dilepton cross sections dcr^^^ by using the EKS98-parametrization of Rf{x, Q"^) ||T2| . 
EKS98 is based on the DGLAP analysis of Ref. [T^ where the nuclear parton distri- 
butions are assumed to be factorizable above Qq = 1.5 GeV. The existence of nuclear 
effects at this scale is taken as a given fact (i.e. the origin is not discussed) and the 
absolute nuclear parton distributions are evolved to the region Q > Qq by using the 
standard lowest order leading twist DGLAP equations. Constraints for Rf{x, Q"^) are 



given by the data on the ratios of structure functions F^^ jF2'^ from deep inelastic 
lepton-nucleus scatterings, by the data on Drell-Yan dilepton production in proton- 
nucleus collisions, and by conservation of baryon number and momentum. More 
detailed discussion of the the analysis 



In p 



is found in [Q. 

it was shown that the ratios i?^ = /^y^/Zf are within a few percent the 
same for different sets of LO parton distributions of the free proton. However, as the 
DGLAP analysis was only performed in the LO, we have a slight inconsistency in 
the NLO computation of the cross sections here. The difference to the NLO DGLAP 
analysis of the nuclear parton distributions is nevertheless expected to be small in 
the ratios -R^^. 



Sensitivity to Parton Distribution Functions 




Figure 2: NLO Drell-Yan rapidity distributions at ^/i 
for various parton distribution functions. 



200 GeV in p-p collisions 



To estimate the sensitivity of the Drell-Yan calculations to the choice of parton 
distribution function sets, we have utilized the most recent analyses of the three major 
groups (MRST99, CTEQ5, and GRV98), where we have used the default parameters 
in the cases that more than one option is available. For further comparison, we 
have also used the previous versions from these same groups (MRSG, CTEQ4, and 
GRV94). 

Shown in Figure ^ are the NLO differential cross sections in center of mass rapidity 
for proton-proton collisions at a/s = 200 GeV at M = 2 and M = 10 GeV. It is 
evident that at M = 10 GeV the x-values probed are sufficiently large such that all 



distribution function sets are constrained by data to very similar values. (Recall that 
at LO the partonic x-values xi and X2 are fixed at xi^2 = On the contrary, at 

M = 2 GeV parton distributions at the lower x- values are not sufficiently constrained 
by data, and variations of 20 - 40 % between sets are not uncommon. Calculations 
with intermediate masses reveal that this large uncertainty essentially disappears 
already when M = 4 GeV, such that above this mass the variations are bounded 
typically by 10 %. 

Corresponding results for ^/s = 5.5 TeV are shown in Figure ^. A similar pattern 
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Figure 3: NLO Drell-Yan rapidity distributions at ^/s = 5.5 TeV in p-p collisions for 
various parton distribution functions. 

exists also at this energy, but with somewhat different magnitudes. At M = 2 GeV, 
the variation is at least a factor of 2. (Note that the apparent discontinuity in slope 
for the CTEQ curves are due to an absolute cutoff below a minimum x-value of 10~^.) 
As one increases M this uncertainly again decreases rapidly, reaching the 10-15 % 
level at M = 10 GeV. 



Shadowing effects in rapidity distributions 

For the effects of shadowing at RHIC and LHC, we choose a representative parton 
distribution set, the MRST99. Shown in Figure ^ are Au-Au LO and NLO cross 
sections per nucleon at RHIC for M = 2 GeV and M = 10 GeV, both with and 
without shadowing. One sees that at smaller M the effects of shadowing are quite 
large, and in much of the kinematic range the magnitude is similar to that of the 
NLO corrections. 




Figure 4: Shadowing in LO and NLO Drell-Yan rapidity distributions calculated with 
MRST99 structure functions for Au-Au colhsions at RHIC with ^/s = 200A GeV. 




Figure 5: Shadowing in LO and NLO Drell-Yan rapidity distributions calculated with 
MRST99 structure functions for Pb-Pb colhsions at LHC with = 5.5A TeV. 
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Figure 6: Ratio of Shadowing/No Shadowing in LO and NLO Drcll-Yan cross sections 
dcr/dMdy as a function of rapidity for fixed values of M at RHIC with = 200A 
GeV. 




Figure 7: Ratio of Shadowing/No Shadowing in LO and NLO Drell-Yan cross sections 
dcr/dMdy as a function of rapidity for fixed values of M at LHC with ^Js — 5.5A 
TeV. 



Figure ^ presents the same calculations at LHC for Pb-Pb collisions. The mag- 
nitudes of the shadowing corrections are even larger here, due to the smaller target 
x-values probed at this energy and rapidity range. 

We also show in Figures ^ and |^ the ratio of cross sections with and without 
shadowing corrections. The dashed lines for LO are determined entirely by the fixed 
Xi^2 -values and the corresponding quark and antiquark shadowing ratios Rf{x, = 
M^) for each nucleus. 

The solid lines for the NLO ratios are not identical to those for LO, but are 
seen to follow those curves quite closely. This behavior is somewhat unexpected, 
since the NLO subprocesses involve integrals over x-values significantly different than 
the Xi^2- It cannot be explained by simple magnitude arguments, since the NLO 
ii'-factors are substantial, and also involve gluon structure functions not present in 
the LO contribution. What is evidently happening is that the NLO integrals can be 
approximated by mean values of the integrand at effective x- values close to the Xi^2- In 
any event, it appears that one may be able to extract overall multiplicative shadowing 
factors which will be approximately independent of the order of perturbative QCD 
and the input structure functions. 

Mass distributions 

We now calculate the expected mass distributions of the NLO nucleus-nucleus 
cross sections at RHIC and LHC, integrated over rapidity intervals appropriate for 
acceptance of the detectors (1.1 - 1.6 at RHIC and 2.0 - 4.0 at LHC). The results 
for RHIC are shown in Figure | for the MRST99, CTEQ5, and GRV98 structure 
functions. Both shadowing and no shadowing calculations are presented. Although 
there is a considerable difference between results with different structure functions, 
especially at low mass, this uncertainty is less than the expected differences between 
shadowing and no shadowing. Also shown are the ratios of shadowing to no shadowing 
predictions for calculations using each of the three structure functions. As anticipated, 
this results in a universal shadowing curve, approximately independent of structure 
function. The shadowing curve at the bottom of the figure actually displays all three 
results - the differences are smaller than the width of the lines. 

Corresponding results for LHC are shown in Figure y. The patterns are very 
similar. At the lowest mass there is somewhat more dispersion due to choice of 
structure function, but we already know that in this region both the reliability of 
pQCD at NLO is quite weak and the parton distributions are less constrained. The 
shadowing ratios converge again to a universal curve. 

For archival purposes, we present in Table |I| the calculated NLO nucleus-nucleus 
cross section values for all three choices of structure functions, both with and with- 
out shadowing. In addition, we list the corresponding p-Au (at RHIC) and p-Pb (at 
LHC) NLO cross sections for the same detector rapidity interval in both the proton 
and nucleus directions. (At ALICE we take into account the asymmetric beam en- 
ergies, and assume that the rapidity interval in the nucleus direction is reached by 
interchanging the beam directions.) 
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Figure 8: Rapidity-Integrated l.l < y < 1.6 NLO Drell-Yan Mass distributions for 
Au-Au collisions at RHIC with = 200A GeV. 
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Figure 9: Rapidity-Integrated 2.0 < y < 4.0 NLO Drell-Yan Mass distributions for 
Pb-Pb collisions at LHC with y/^ = 5.5A TeV. 



Conclusions 



We have studied tlie production of Drell-Yan dileptons at energies and phase space 
regions appropriate for the BNL-RHIC and CERN-LHC heavy ion experiments. The 
NNLO calculations, available for da/dM, indicate that a perturbative calculation 
converges well for M > 4 GeV, the NNLO correction being < 5 %. At the LHC 
energy it becomes rapidly less reliable as M approaches 2 GeV. Using a representa- 
tive selection of most recent parton distribution sets, MRS991 CTEQ5 0, and 



GRV98 |Tn], we found that the differences between the sets lead to an uncertainty 
in the cross sections which at small masses, M ~ 2 GeV, is typically 20 - 40 % at 
RHIC but around a factor of two at LHC. This reflects the lack of precision exper- 
imental constraints on parton distributions at small x. At RHIC the results from 
different sets converge as M increases, but at LHC persist at the 10 - 15 % level even 
for M = 10 GeV. We investigated the nuclear effects using the parton distribution 
functions for free nuclei modified with multiplicative shadowing functions determined 
from the available experimental information [jl2[. At RHIC energy, shadowing re- 
duces the cross section at M ~ 2 GeV by a universal (i.e. independent of structure 
function set or order or perturbative QCD) factor of ~ 35 %. This factor approaches 
~ 10 % for M ~ 10 GeV. At LHC the reduction for small masses is ~ 55 % and 
more than 30 % even at M = 10 GeV. In summary, these calculations are limited by 
uncertainties due to convergence of the perturbative series and parton distributions 
which become large at small values of M/ y/s. Thus they are generally not a problem 
for all of the RHIC results, but become severe at LHC for M < 4 GeV. The overall 
effects of shadowing at RHIC energies are smaller in magnitude than those at LHC, 
but the change of the mass distribution shape due to shadowing is a larger effect at 
RHIC. 
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Table 1: NLO Drell-Yan Calculations for RHIC and LHC 



